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ABSTRACT: Virtual fermionic Nf = 1 and Nf = 2 contributions to Bhabha scattering are combined 
with realistic real corrections at next-to-next-to-leading order in QED. The virtual corrections are 
determined by the package BHA_NNLO_HF, and real corrections with the Monte Carlo generators 
Bhagen-IPh, Helac-Phegas and Ekhara. Numerical results are discussed at the energies 
of and with realistic cuts used at the <I> factory DA<I>NE, at the B factories PEP-II and KEK, and at 
the charm/t factory BEPC II. We compare these complete calculations with the approximate ones 
realized in the generator BabaYaga@NLO used at meson factories to evaluate their luminosities. 
For realistic reference event selections we find agreement for the NNLO leptonic and hadronic 
corrections within 0.07% or better and conclude that they are well accounted for in the generator 
by comparison with the present experimental accuracy. 



KEYWORDS: Electromagnetic Processes and Properties, NLO Computations. 



Contents 



1. Introduction 1 

2. The NNLO massive corrections 3 

2. 1 The virtual plus soft photon corrections <J^° 5 

2.2 Hard photonic corrections a^ NL0 9 

2.3 Real electron pair contributions a g L ° g _^ +e _^ 11 

2.4 Real muon and tau pair contributions o^+^+^n 12 

2.5 Real pion pair contributions C^ e _, +J[ _^ 13 

2.6 Other hadronic corrections 14 

2.7 Exact NNLO numerical results 17 

3. The NNLO massive corrections in B abaYaga @ NLO 19 

3.1 The program 19 

3.2 General formulation of Bab aYag A @ NLO 20 

3.3 NNLO massive corrections in BAB AYAGA@ NLO 22 

3.4 Numerical results of Bab aYag A @ NLO 24 

4. BabaYaga@ NLO versus the exact NNLO massive corrections 25 

5. Conclusions 26 

A. The experimental cuts 32 

A. 1 <I> factory KLOE/DA<I>NE (Frascati) - the reference event selection 32 

A.2 charm/t factory BES-III (Beijing) - the reference event selection 32 

A.3 B factory BaBar/PEP-II (SLAC) - the reference event selection 32 

A.4 B factory Belle (KEK) - the reference event selection 33 

References 34 



1. Introduction 

The Bhabha scattering process 

e + e -^ e + e - (1.1) 

is an invaluable tool for the luminosity determination in various experiments. Both low energy 
devices, operating from about 1 GeV to several GeV and high energy devices, planned to operate at 
hundreds or thousands of GeV, require theoretical predictions for the Bhabha cross section with quite 
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accurate determinations of QED radiative corrections. The latter contain, besides exponentiated 
leading logarithmic terms, also the complete fixed order contributions, and in particular the complete 
two-loop QED corrections. 

Aiming at per mille accuracy or slightly better, the radiative corrections may neglect the constant 
terms in the electron mass m e . At next-to-leading order (NLO) final states with unresolved photons 
will contribute, 

e + e~ -> e + e~(y). (1.2) 

Further, new mass scales start to play a role, but only in the one-loop self-energy insertions. At 
next-to-next-to-leading order (NNLO), further final states of Bhabha scattering contain unresolved 
photons, fermion pairs, or hadrons: 

e + e~ — > e + e~(y,YY), e + e~(e + e~), <? + <?~(/ + /~), e + e~ (hadrons) . (1.3) 

At this order of perturbation theory, a variety of Feynman diagrams depend on additional mass 
parameters, and one may formally distinguish between Nf = 1 corrections (with only electrons) and 
Nf = 2 corrections, and the latter ones are technically more complicated due to the additional mass 
scale. 

The NLO corrections, with inclusion of certain leading higher order terms, are known since a 
while and several Monte Carlo (MC) programs are carefully tuned. A recent comprehensive review 
on precision predictions for scattering experiments at meson factories contains a detailed discussion 
of the state of the art [1]. 

As far as virtual corrections are concerned, in the last few years there has been major progress 
in the evaluation of the corrections at the NNLO accuracy. In fact, the photonic two-loop QED 
corrections were first evaluated in the massless case in [2] . The photonic corrections to massive 
Bhabha scattering with enhancing powers of \n(s/m 2 ) were soon derived from that [3]. The missing 
constant term in m e [4] plus the corrections with electron loop insertions [5, 6, 7, 8], called the 
Nf = 1 case, followed not much later. The heavy fermion (or Nf = 2) corrections were first derived 
in the limit m 2 « m 2 << s, \t\, \u\ [8, 9], where nif is the mass of the heavy fermion and s,t,u are 
the usual Mandelstam variables, and soon after also for m 2 « m 2 ,s, \t\, \u\ [10, 11]. Finally, using 
dispersion relations, also hadronic corrections became known [12, 13, 14]. 

A complete collection of all the relevant formulae for the massive virtual NNLO corrections 
used in this paper can be found in the just mentioned papers and in [15]. 

When fermion loops or virtual hadronic corrections are taken into account, the question of 
considering also the real emission of the corresponding particles arises. That was studied for 
the emission of electron pairs in [16] in the soft limit of electron pair energy and in logarithmic 
accuracy. It is shown that the leading logarithmic corrections In 3 (s/m 2 ) cancel with those from the 
irreducible two-loop vertex corrections with electron loops. A similar cancellation is expected for 
the combination of heavy fermion pair emission with irreducible two-loop vertex corrections with a 
heavy fermion loop. In practice, however, the situation is evidently a bit more involved, especially 
at smaller energies, when s, \t\, \u\ ~ m 2 . Then, the logarithms are not numerically dominating and 
more diagrams get important. Nowadays, MC programs can do that job. In this article, we will 
perform such a study of reaction (1.3) due to the additional emission of real pairs of leptons with the 
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Fortran packages HELAC-PHEGAS [17, 18, 19, 20], similarly to what was done in the 1990s for 
small-angle Bhabha scattering at LEP [21, 22]. For heavy fermions and hadrons, we present here 
the corresponding results for the first time. The case of real hadron emission in Bhabha scattering 
deserves special attention. The only existing event generator contains only pion pair emission and 
the results obtained for the real pion pair emission serve as an indication of the size of other left 
over corrections. For a consistent treatment, we replace the hadronic dispersion integrals in the 
virtual and soft real hadronic corrections as described in [12, 13] by the pion pair form factor. The 
prediction is then combined consistently with real pion pair emission as evaluated with the MC 
event generator Ekhara [23, 24, 25, 26]. In the calculations we use the pion form factor from [27]. 
For the virtual and hard photon corrections the full hadronic corrections were obtained using the 
vacuum polarisation insertions. 

The results may be compared with those from the Bhabha generators which are usually applied 
for experimental simulations of Bhabha scattering; here we look at BabaYaga [28, 29, 30], in 
particular to the latest and most accurate version BabaYaga@NLO [31]. 

The aim of the article is to put together all the above discussed NNLO corrections to Bhabha 
scattering taking into account real experimental conditions and examine how well they are accounted 
for in the event generator BabaYaga@NLO used at meson factories for their luminosity measure- 
ments. So far, at NNLO level, virtual corrections have been checked in detail only for situations 
where the dependence of soft radiation on the maximum soft photon energy ft) (or, equivalently, 
the minimal hard photon energy) is "switched off" by setting ft) = yfs/2 [4, 12, 13, 14]. This was a 
good way to compare results obtained by different theoretical groups, but certainly has nothing to 
do with reality. We restrict ourselves here to low energies (meson factories) because presently they 
are of immediate relevance from the experimental point of view. 

We just mention for completeness the last remaining NNLO issue: that of radiative loop 
corrections, i.e. the NNLO contributions from the interference of photonic bremsstrahlung off 
one-loop diagrams with lowest order real photon contributions, first studied in [32] with a restriction 
to the factorising diagrams. The technical complications arise from non-factorising diagrams, the so- 
called pentagon diagrams. Recent papers on this issue are [33, 34, 35, 36], but so far without explicit 
numerical results. Sample numbers for the virtual one-loop (plus real soft) QED corrections to the 
hard-bremsstrahlung process e + e~ — > e + e~yare given in [36]. For future measurements, it would be 
worthwhile to answer the question if (and when) these corrections have to be yet included in the MC 
event generators employed for simulating Bhabha scattering events at low-energy high-luminosity 
electron-positron colliders. 

The paper is organised as follows. In Section 2 we discuss the exact NNLO massive cor- 
rections to Bhabha scattering and present benchmark results for event selections close to the 
experimental ones. In Section 3 we describe the approximate treatment of these corrections in the 
BabaYaga @NLO event generator and derive benchmark results for the same event selections 
as for the exact results. In Section 4 we show detailed numerical studies of the quality of the 
approximations used in BabaYaga @NLO. We draw our conclusions in Section 5. 

2. The NNLO massive corrections 

The complete NNLO Nf = 1,2 corrections to Bhabha scattering consist of three parts, each of them 
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(a) (b) (c) (d) 

Figure 2.1: (a)—(c) are sample two-loop diagrams; their interference with (d) is contributing to 
a^ LO , part ofa^°, Eq. (2.1). 



with contributions from virtual and real electron pair corrections (Nf = 1 case) and corrections due 
to muon pairs, T pairs and hadrons (Nf = 2 cases): 

da da da da 

_ do e+e - do^+n- da T + T - do had 

~ da da da da ' { 

We want to concentrate here on the interplay of virtual and real corrections. For the various pure 
self-energy corrections in CJ l ™ LO we refer to [1, 13] and the references quoted therein. This is in 
accordance with the approach chosen in the MC packages used for the interpretation of experimental 
results, and we will not include these pure two-loop self-energy corrections in the numerical results 
discussed below. 

As a result, the following contributions will be studied: 

• the CJ^ L0 consists of virtual two-loop corrections shown in Fig. 2.1 and loop-by-loop 
corrections <t]5jl shown in Fig. 2.2: 

^.NNLO ^.NNLO , ,-NNLO o\ 

°virt — °2L T °7ZJL \ L - L > 

• contributions with real photon emission, shown in Fig. 2.3: 

• contributions with real pair or hadron emission depend a bit more on the flavour, as shown in 
Figs. 2.4-2.6: 

°real = °e+e- (e+e~ ) + °e+ e - (/+/- ) + °e+ e - {hadrons) (2.4) 

The self-energy blobs in Figs. 2.1-2.3 stand for lepton pair or hadronic self-energy insertions. 
Technically, we have to re-order the corrections. The following pieces will be summed up: 
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(a) (b) (c) 

Figure 2.2: (a)-(b) are interfering with(c); they are samples of the so-called loop-by-loop corrections 
a™ L ° part qfo%%°, Eq. (2.1). 

• virtual plus soft photonic corrections: 

oZ° = <T + ^ ft ((0), (2.5) 

see Figs. 2.2-2.3. The sum is infrared finite, but depends on a soft-photon cut-off parameter 
co; 

• hard photon radiation: 

„NNLO _ . 

'y.hard 



NNLO— NLO (CO). (2.6) 



Here we take into account realistic experimental phase space cuts. The sum cr^ + cr^ NL0 
now is independent of CO. 

• Bhabha scattering with the additional production of fermionic pairs or of hadrons, o^° al . 
2.1 The virtual plus soft photon corrections <J^° 

The Feynman diagrams of the virtual corrections contain electron, heavy lepton or hadronic self- 
energy insertions. We include in the cross section (i) the irreducible vertex corrections of Fig. 2.1(a) 
(see Eq. 39 in [13]); (ii) the factorizable vertex and box corrections of Fig. 2.1(b) and Fig. 2.2(a,b) 
(see Eqs. 58, 60, 61 in [13]); (Hi) the irreducible box corrections of Fig. 2.1(c) (see Eq. 65 in [13]); 
(iv) the soft real photonic corrections (not shown, see Eqs. C4 and 62 in [13]). For electrons, we 




(a) (b) 
Figure 2.3: Interference of (a) and (b) is a sample contribution to Oy L0 in Eq. (2.1). 



will take the exact expressions, and for the other contributions the approximation m 2 e « s,\t\,\u\ is 
applied. 

The irreducible vertex diagrams are infrared finite, but the reducible vertices and the box 
contributions are not. To make the latter two infrared finite, one has to add the corresponding soft 
real photon emission: 

da%? doS£ r da^ f - doSgj do%? 

~doT = do. + / ^ T da + + x F r**W (2 " 7) 

Here, the term a^f x Fy. so f t ((o) arises from the interference of (soft) single-photon bremsstrahlung 
diagrams, where one of the diagrams has a self-energy insertion. The a,^° comes also from a 
Bhabha cross section due to diagrams with first-order fermionic (hadronic) self-energy insertions, 
and Fy^ f t (co) is the usual soft photon eikonal factor; for explicit expressions see [8, 13]. The 
parameter ft) is the infrared cut-off 

a = E yjqfi = E rhard ( 2 - 8 ) 

and has to be adapted such that soft-photon emission has the Born kinematics and the sum of soft 
and hard photon radiation is numerically independent of ft); typically, it is co/E\, eam = 10~ 6 • • • 10~ 3 . 
The evaluation of the NNLO virtual corrections has been detailed elsewhere and we may restrict 
ourselves here to few remarks on the specifics of this article. 
In the simplest case of a one-loop self-energy insertion, 

^n (s »)(,„-^), (2.,, 

one may just evaluate Feynman diagrams and gets after renormalisation for the case of a light 
fermion with mass m,m 2 << s: 



5 , / m 
+ ln 



,2 



3 \ s + iS 



(2.10) 



For the virtual hadronic and heavy lepton pair corrections we use the dispersion approach. Here, the 
photon propagator is substituted in the Feynman diagrams as follows: 

guv irff 2\ { a rdzR(z) ,2\{ 



q 2 + i5 q 2 + i8 J 3k J m2q z -W^jy^v q2 + i8 

with the propagator 

K SE (q\z)= 2 1 ■ (2.12) 
q l - z. + 10 

For heavy leptons with mass m/ and charge Qi = —1, it is at one loop (exact in the mass m/): 



l + 2-MWl-4-^. (2.13) 
For hadronic corrections, a natural choice for RhM is a representation by experimental data: 



-6- 



where Ghad(z) = c(e + e~ — > y* — > hadrons;z). The real hadronic emission can be studied at present 
only for pion pair production as only for this hadronic final state a MC code is available. Corre- 
spondingly, for the pion case we use instead of (2. 14) the undressed pion form factor F K \ 



Ohad(z) -> a(e + e- -> n+n~) = ^°^l\F K (z)\ 2 , (2.15) 

3 Z 



where 



p n = (\-Am 2 Jz) xl2 , (2.16) 
is the pion velocity. For the cross section ratio R, which we need here, this transforms to: 



Rhad(z) -> Rxn{z) = fiFM)] 2 . (2.17) 



The pion form factor F n (z) has been determined in [27]. The numerical studies using this parameter- 
isation are presented in Section 4. 

All this applies to one-loop insertions in reducible diagrams. For irreducible two-loop vertex 
and box diagrams, one has to perform an additional loop integration, and the result is more involved. 
The complete virtual NNLO Nf = 1 corrections (due to electron self-energy corrections) 0^°+ e - 
are known exact in m e [7] in terms of harmonic polylogarithms [37]; a re-calculation and also 
corresponding formulae in the kinematic limit m 2 . « s,t, in terms of ordinary polylogarithms, are 
given in [8]. The virtual NNLO Nf = 2 corrections (due to hadronic self-energy corrections) 
are determined with dispersion formulas. 

At the end of this Section, we shortly comment on the structure of the various contributions to 
the virtual plus soft photon cross section: 

NNLO _ NNLO, NLO / \ 
°v+s — a virt + O y,soft\ w ) 

_ rt NNLO j_ „NNLO , „NNLO , ^NLO v 77 / ,,,\ n ,m 

- °fact +°vert + °box + G virt X 1 Y,soft\G>) ■ (2.1%) 

The eikonal factor is in the limit of small m e [13]: 



a 

Fy,soft(G>) = — 



ln(V<)-21n v _ 



21n (s/m 2 e ) +2-21n ' 
In (s/m 2 e ) 2 -4^2 + 2 In (s/m 2 e ) +2Li 2 (~) - 2Li 2 (-") | 



(2.19) 
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The virtual contributions are: 



_NLO 



J fact 



_NNLO 
°vert 



J box 



a\3 



71, 

a\ 4 

71 

a 



~ ( - I < o, 

71 



+ *E L^^[C L,, ^(-v-<-)+C L 'Xs£(>..-)] ^. (2.20) 



M ■'^0 



J fact,e+e- 



yNNLO 

vert,e+e" 



*E/ 2 



<1zRm{z) 



[C^K S E(s,z)+C^°K SE (t,z)] , (2.21) 



+*EL 



M J «0 



<fe jz) r„ NN LO y (c I /""NNLO y ft A] 
Invert, s Invert \ A iZ) <^vert,t ^vert \ l ; Z)\ 



-NNLO 

box,e+e~ 



dzR-M{z) 



{Cboxs [ K box,A{s,t,z)+K hox £(t,S,z) 



+ K box ^ c (u,t,z) -K boxM (u,s,z)] +C™ x ° t [K hoXtB (s,t,z) +K boxA (t,s,z) 



K h0X:B (u,t,z) + K h0XjC (u,s,z)} | 



(2.22) 



(2.23) 



The sum over M covers m ll ,m x ,m n with the corresponding parameterisations of Rm(z); for leptons 
see (2.13). The lower integration bound is Mq = Am 2 for leptons. For hadrons, where one sums 
over all the hadronic contributions, the lower bound is M 1 = m^, corresponding to 7T°y, the lightest 
hadronic final state. 

The kinematical factors C (rational functions of s and t) and the kernel functions K are universal. 
A special role play the irreducible vertex and box diagrams. In these diagrams, the self-energy 
correction is part of a loop insertion, and due to the additional loop momentum integration the 
replacement (2.1 1)-(2. 12) leads to more involved kernel functions compared to (2.12); for the vertex 
[38]: 



K vert (x;z) 



I/_Z 

3 1 8 



x 

z. 



Here hi2(x) is the usual dilogarithm and £2 = Li2(l) = 7C 2 /6. From the irreducible box diagrams 
we have three different, lengthy box kernel functions K box ^(x,y, z) , K box ^(x,y, z) , K boXt c(x,y, z)', see 
for explicit expressions Eqs. (71)-(73) of [13]. 

For practical reasons, it makes sense to split the <J^j° into two pieces, namely the infrared 
finite cO° of (2.22) and the so-called "rest" a r N ^°, 



—NNLO —NNLO , —NNLO 



(2-25) 



We just remind the reader that a third piece, the pure self-energy corrections, is not included in 
the study. The cr, N e ^° is the sum of all the infrared divergent contributions. In [13], it is detailed in 
Section VI and in Eq. (87). This sum is infrared finite, but depends on the photon energy cut co 
related to the separation of soft and hard photons. In the energy regions of relevance with s > Mq, 
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the net result may adapted from Eq. (93) of [13]: 



7l 2 S 



*«W-i-«<«) 

Mk Z t-z 



(2.26) 



+ rdz^-^{R(z)F 2 (z)-R(s)F 2 ( S ) + [R(z)F,(z)-R(s)F 3 { S )]\n 1-M 



+ 



R(s) 



{f 2 (s) ln( 



M 2 



+ F 3 (s) 2£ 2 + - In 



M 2 



■6C 2 F 4 (j) 



1 +Li 2 1 



Ml 



The explicit expressions for F\ to F4 are given in (88)-(91) of [13]. They are infrared finite, but 
depend on 2(0 f^Js. 

It is well-known that the irreducible vertex corrections from a fermion pair with mass m con- 
tribute to the cross section with terms of order In 3 (s/m 1 ). For electrons, this is a huge enhancement: 



and the form factor is for m 2 e /x « 1 [39] : 



V 2 e(x): 



1 

36 



hr 



H In" 

72 



+ v 2e 



(2.27) 



(2.28) 



The vertex function may be found in [40] exact in m e . For heavy leptons, the logarithmic terms in 
V 2 f(x) agree with V 2e (x), but a deviation appears in the constant term which becomes v 2 f [39, 8]: 



Vie 



V 2 f 



+ C2 In 



+ 



1 /265 
6 

1 /3355 19 „ 



1 /383 
27 



C2 



(2.29) 
(2.30) 



It is these logaiithmic dependences which make the real pair emission contributions so important, 
because cancellations of the leading terms appear. 

Technically, we evaluate the electron corrections dc™\° e /dQ. with the Mathematica program 

, NNLO v+s c 

CROSSSECTION.M [15]. The contribution aa da is represented there by function NNLOel. 
However, it includes also the iterated one-loop self-energies NNLOFE2 and the genuine two-loop 
self-energy NNLOfeI. Both are calculated separately, but are not included here in the numerics. 
They have been subtracted from NNLOel in order to estimate what we call here cr^°. The heavy 
fermion corrections have been calculated with a Fortran package applying the dispersion technique 
described here. In order to cover also pion pair corrections (7,™^, the c^ff is determined with an 
updated version BHA_NNLO_HF of the Fortran package BHBHNNLOHF [13, 15]. 

2.2 Hard photonic corrections a^ NLO 

The NNLO hard photonic corrections a^ NL0 with a self-energy insertion arise from the classes of 
diagrams shown in Fig. 2.3, with emission of one hard photon. They were calculated with the 
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Fortran program Bhagen- IPh-VAC [41] based on the generator Bhagen- IPh [42]. This cross 
section depends on the soft photon cut-off £™ m = ft) and only after adding them to cr^J: , the sum 
of the two h is independent of the cut-off. We calculate here separately the contributions from 
diagrams with the electron, muon, tau, pion and complete hadronic vacuum polarisation insertions. 
The dependence on additional cuts is crucial and varies considerably with the experimental set-up. 
A careful discussion is given in Section 4. 

Even if the insertion of the vacuum polarisation corrections to the square of the tree level 
amplitude [43] is straightforward, for completeness we give below the formulae, which are used in 
the unpublished program [41]. We consider the process 



e + {p+)+e (p-)^e^ 



+ e (q-)+y(k) 



(2.31) 



and follow the notation of [43] 



s= (P++P-) 2 , t = (p+-q + ) 2 , u= {p+-q-) 2 , si = {q + +q-) 2 , 
h = (p- - q-) 2 , «i = (p- - q+) 2 , k±=p±.k, h±=q±.k. 

The differential cross section for the process (2.31) can be written as 



do 



a 3 , d 3 q+d 3 q-d 3 k 4 

^(X + F+Z) — — — 5 (p ++P --q + -q--k), 



(2.32) 



(2.33) 



where E + ,E-, E Y are the energies of the final positron, electron and photon, respectively. The 
quantities X, Y, Z refer to the ^-channel annihilation, the ^-channel scattering and the interference 
part of the squared amplitude, respectively. Keeping only the diagrams with a virtual photon 
exchange (the original formulae [43] contained also the weak Z boson contributions) and with the 
vacuum polarisation corrections included they read: 



X = (Re(U(s)+n(s l )))(t 2 + tj + u 2 + ul 
+ 2R<z(n{s 1 ))({t 2 + t 2 + u 2 + u 2 ^ 



1 



4s si 



+ 



k + h- k-h + k + h- 



k_h_ 



2Re(n( J ))^§ 



■ + 



4s\k + k 



+ ■ 



+ 2Re(U(s))({t 2 + t 2 + u 2 + uj) 



4sh+h- 



+ 



(h + ) 2 (h.) 2 (h + ) 2 (/*_) 



2Re(n(si)) 



2s\ 



+ 



+ 



■ + ■ 



{k + f {k^f (k + ) 2 



(2.34) 



Y = [(U{t) + U{t l )){s 2 + s 2 l +u 2 + uj)] 



4th 



U U] S S] 



k + h- k-h 



k + k- h + h- 



2Yl{t){s 2 + s\ + u 2 + u\ 



^-m^w^+u 2 )^- 



-2U{h) 



2t\ 



+ 



+ 



■ + ■ 



(h + y {k + f (h + ) 2 (k + ) 



2U(t) 



2t 2 



+ 



+ ■ 



+ ■ 



(/i_) 2 (k_) 2 (h^) 2 (jfc_) 



(2.35) 
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(a) 



(b) 



(c) 



Figure 2.4: Samples of the 36 diagrams contributing to e + e — > e + e e + e 



z u 2 + u\ [n(/)+Re(n(j)) ( u 



St 



+ 



k-h + h+-h- k-h- 



| n(f!)+Re(n(g)) / ui | s 



st I 



t\ 



kih- h + h- k+h.4 



| n(Q+Re(n(^i)) / u y | si 



sit 



k+h- k+k- k-h- 



n(fi)+Re(n(ji)) f u si t\ 



S\t\ 



+ 



k-h+ k+k- k+h 



2 

171 U 

-Mn( fl )+Re(n(,)))-^ 

st i (h+) 



+ 

2 2 

' ^O+Re^)))^ 



(*+) 2 



- £ m) + Rc(n W) ) ^ - =f W ) +Re(n (jl ))) ^ • 



(2.36) 



The generation of the phase space was not changed with respect to the original program BHAGEN- 
IPh [42]. 



2.3 Real electron pair contributions a L ° + _s 

The most important real fermion pair corrections to Bhabha scattering are, at all energies, the 
unresolved electron pair corrections. There are 36 diagrams of this kind contributing to ^+ e -( e + e -y 
part of the Bhabha cross section (2.1). Sample diagrams are shown in Fig. 2.4. 

The e + e~ pair corrections fall into three classes: four ^-channel diagrams with two e + e~ 
pairs in the final state, 24 diagrams (8 in ^-channel and 16 in ^-channel) with one e + e~ pair, and 
8 peripheral ^-channel diagrams (no e + e~ pair). What is usually considered as the electron pair 
corrections, are those with two electron pairs in the final state, Fig. 2.4(a). The contribution from 
such soft electron pairs is known [16], see also [13]. It is, in the limit of small m e and small energy 
cut-off parameter D of the unresolved e + e~ pair, proportional to the lowest order Bhabha Born cross 
section a4? _ : 



-{e+e-) ~ °e+e- 



(2.37) 
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with: 



"soft 



l -L] +L 2 S (2ln(D) - ^ +L S U\n 2 (D) - y ln(D) +A, j (2.38) 



where 



+ l -I? t +L 2 f 21n(D) - +Lt ^41n 2 (D) - y ln(D) +A, 
- hi -Ll (l\n{D) - —L u Uln 2 (D) - y ln(D) +A„ 



L s = In ( -U , (2.39) 

L v = In ( — % ) , v = t,u, (2.40) 

V m eJ 

A, = y-4C 2 , (2.41) 

/l±COS0\ 

A v = A,. + 2Li 2 , v = t,u. (2.42) 



The parameter D has to fulfill: 



2m e < < DE beam < < E beam . (2.43) 



From the sum of (2.37) and (2.27), the compensation of the leading mass singularities (contained 
here in the L^,L ( 3 ,L^ terms) in the cross section becomes evident. 

If there are unresolved contributions, which do not fulfill Born kinematics, or if the logarithms 
are not really big, or if the experimental accuracy is at the per mille level or better, a complete 
calculation is needed, and this is part of the present study. 

The Feynman diagrams are finite as long as the electron mass is assumed to be finite, so that a 
straightforward Feynman diagram calculation of the 2 — > 4 process can be performed without any 
true singularities. 

2.4 Real muon and tau pair contributions <7g+ e (/+/ j 

For both e + e~ — > e + e~li + [l~ and e + e~ — > e + e~T + T~ there are 12 diagrams. Samples of them are 
shown in Fig. 2.5. 

There are four classes of diagrams, to be discussed here for muon corrections: two ^-channel 
diagrams with production of both an e + e~- and a fi + pL -pair; two ^-channel diagrams with an 
e + e~-pair; six diagrams (two in ^-channel, four in ^-channel) with a fermion pair; and finally two 
peripheral r-channel diagrams. Contrary to real electron pair corrections, for heavy lepton pairs it 
is at the meson factory energies never appropriate to assume that the lepton mass is much smaller 
than e.g. yfs. So, one has to perform a complete lowest order 2 — > 4 Feynman diagram calculation. 
Again, the cross section is finite as long as all masses are assumed to be finite. The same discussion 
holds for the process e + e~ — > e + e~r + T~. 

The results for the electron, muon and tau pair corrections to the Bhabha scattering process 
have been obtained in the framework of the HELAC-PHEGAS leading-order MC program [17, 20]. 
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(a) (b) (c) (d) 

Figure 2.5: Samples of the 12 diagrams contributing to e + e~ — > e + e~}X + }X~. A similar set of 
diagrams describes e + e~ — > e + e~x + T~. 

The phase space integration was executed with the help of PHEGAS [18], a general purpose multi- 
channel phase space generator. Helac-Phegas generates, in a fully automatic manner, events for 
all possible parton level processes at hadron and lepton colliders within the Standard Model. More 
precisely, integrated cross sections and kinematic distributions with arbitrary cuts on particles in the 
final state and with full spin correlations can be obtained. It has already been extensively used and 
tested in phenomenological studies, see e.g. [19, 44, 45, 46]. 

In the present study, the exact QED 2—^4 matrix elements for e + e~ — > e + e~£ + £~ processes, 
where ^ = e ± ,/j. ± ,z ± , have been generated, including all Feynman diagrams (36, 12 and 12 
respectively) and mass terms. Let us mention here, that in order to generate pure QED contributions 
in Helac-Phegas the coupling of £ + £~ to the Z boson has to be simply set to zero. We have 
checked that the Z contributions are negligible for event selections used in this paper and do not 
affect any conclusions. 

2.5 Real pion pair contributions a™ e -, n+n -s 

The lightest hadronic final state produced in e + e~ scattering via the one photon exchange mechanism 
is the charged pion pair. This final state (i.e. e + e~ n + K~) was investigated in details in [24, 25, 26] 
and implemented into the Monte Carlo generator E KHAR A [24, 23]. Since other hadronic final 
states produced this way were never implemented in a generator, their studies are impossible at 
present and results obtained for the reaction e + e~ — > e + e~7t + 7t~ will be used also as a hint towards 
understanding the importance of other reactions like e + e~ — > e + e~ + (7r + 7r~7r°, K + K , KsKl, ■■■). 
The corresponding set of diagrams consists of 14 diagrams, with their representatives shown in 
Fig. 2.6. To model the pion-photon interactions we use the vector dominance model with the pion 
form factor from [27]. These contributions can be seen as: initial state electron pair emission (a), 
final state electron/pion pair emission (b,c), pion pair emission from the ^-channel Bhabha process 
(d) and y* — 7* pion pair production (e). The last set of diagrams is small for large electron and 
positron angles and its modelling is relatively crude, neglecting scalar meson production and the 
subsequent decays to pion pairs. 
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(a) 



(b) 



(c) 





(d) 



(c) 



Figure 2.6: Sample diagrams with real pion pair emission. 



2.6 Other hadronic corrections 

As discussed in the previous Section the lack of an event generator for processes e + e~ — > e + e~ plus 
hadrons does not allow the computation of the real hadron emission with the exception of charged 
pion pairs. However contributions coming from virtual hadronic vacuum polarisation insertions can 
be calculated completely (see Section 2.7). In this Section we show how the vacuum polarisation 
insertions from pions compare to the full hadronic corrections. It has to be stressed that even if this 
difference can give some indication of the size of the missing real emission contributions, its actual 
size depends heavily on the event selection used by a given experiment. Thus a reliable estimation 
of the missing terms is not possible without Monte Carlo simulations. 

Specific problems are caused by narrow resonances like J/y, y(2S), ... and one should devote 
them a special treatment. Narrow resonances with mass M res and partial width T^/ can be described 
approximately by the ansatz 



Rres{z) 



9k 
a- 



M res Y e re f8{z-M 2 res ). 



(2.44) 



Based on this, their contributions to the NNLO Bhabha process can be derived from the general 
formulae of [13]. We discuss here as an example the contribution from the "rest" (Eq. 2.26); 
according to Eq. (87) of [13] it reads: 



9a 2 r? 



dQ. 



% s M„ 



-M 2 



+ 



1 



-M 2 



F 2 (M r 2 s )+F 3 (M r 2 es )ln 



1 



M 



(2.45) 
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This is basically also Eq. (E4) of [13, 15]. Eq. 2.45 becomes invalid when the center of mass energy 
comes too close to the position of a resonance, i.e. if (s — M r 2 es ) < r^ s e M res . In the numerical 
examples, Table 2. 1 , this is not the case. 

Table 2.1: Soft+virtual NNLO contributions o™\^ es from narrow resonances (n.r.) defined by 
Eq. (2.45) for the Bhabha prOC€SS with CO /-Ebeam — 

l(T 4 (innb). The narrow resonance located 
closest to the center of mass energy of the given collider is included (first column, res) and excluded 
( second column, res'). The third column contains the Born cross section. 





—NNLO 
°rest,res 


—NNLO 
rest.res' 


Ob 


KLOE 1.020 


[all n.r.] 
-0.04538 


[n.r. without J/y/{lS)] 
-0.0096 


529.5 


BES 3.097 


[all n.r.] 
228.08 


[n.r. without J/y/{lS)] 
-0.0258 


14.75 


BES 3.650 


[all n.r.] 
-0.1907 


[n.r. without y(2S)] 
-0.023668 


123.94 


BES 3.686 


[all n.r.] 
-62.537 


[n.r. without y/(2S)] 
-0.0254 


121.53 


BaBar 10.56 


[all n.r.] 
-0.0163 


[n.r. without T(4S)] 
-0.01438 


6.744 


Belle 10.58 


[all n.r.] 
0.04393 


[n.r. without T(4S)] 
-0.0137 


6.331 



To illustrate the role of narrow resonances, in Table 2. 1 we show numerical results based on 
Eq. (2.45). We use parameters listed in Table 2.2. We can see that the contributions from narrow 
resonances dominate the NNLO Bhabha correction for BES running at J/y and y(2S) energies. 
For the remaining cases narrow resonances contribute below the per mille level when compared to 
the Born cross section Gb or to BabaYaga@NLO best predictions a BY , see Table 4.1. 

We conclude that for experiments performed on top of a narrow resonance, this resonance 
cannot be treated as a mere correction and more detailed studies have to be performed. These should 
include examining of finite width effects, beam spread effects, estimation of NNNLO corrections and 
the accuracy of the vacuum polarisation insertions in a close vicinity of these resonances. Having 
this in mind we do not present here hadronic contributions for the BES-III experiment running at 
// \j/ and y(2S) energies and we plan to devote to this issue a separate study. 

We now come to the net hadronic vacuum polarization effects, i.e. look now at the sum of 
the so-called "rest" terms and the irreducible vertex corrections. To obtain all numerical results 
below we use "rest" as given by Eq. 2.26 together with the corresponding formula for the vertex, Eq. 
(2.22). The recent update of R] KU \ valid in the range mi < s < (100 GeV) 2 [47] is applied, and for 
higher s the Rj uu i is taken from [48]. For further details on the implementation of R] w d see Appendix 
Eof [13]. 

The numerical integrations for the hadronic virtual and soft contributions were performed by 
means of the adaptive integration routine VEGAS [49] , which works efficiently even for so narrow 
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Figure 2.7: The R] xa d vs. R n +n -■ Narrow resonances are not included. 
resonances like J/y. 

Table 2.2: Parameters of narrow resonances used in [47] (T. Teubner, private information). 



resonance 


M res [GeV] 


Tf e f [keV] 


Vv(is) 


3.096916 


5.55 


V(2S) 


3.686093 


2.33 


T(1S) 


9.46030 


1.34 


T(2S) 


10.02326 


0.612 


T(3S) 


10.3552 


0.443 


T(4S) 


10.5794 


0.272 


T(5S) 


10.865 


0.31 


T(6S) 


11.019 


0.13 



In Fig. 2.7 we compare the full result for R] KU i and the contributions coming from pions only. At 
low energies (Fig. 2.7 a) the biggest difference comes mostly from contributions of three pions and 
from kaon pairs pronounced at ft) and resonances, while at high energies the pion contributions 
vanish rapidly and do not play any significant role (Fig. 2.7 b). 

From Table 2.3 it is clear that the pion contributions are not sufficient for an accurate evaluation 
of the hadronic contributions and the lack of a generator for generic hadronic final states does not 
allow to draw a final conclusion about the real hadronic corrections with the exception of the charged 
pion pair emission and event selections which kinematically exclude the hadron production (KLOE 
energy). 

The R enters the results with weight functions, which give more relevance to the low energy 
range, however the differences are important also there. Partially one can see the effect of the 
weight functions comparing the full hadronic corrections to the vacuum polarisation with pion pair 
contributions (Fig. 2.8), however the complete weight functions are complicated and different for 
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Table 2.3: Comparison ofhadronic contributions modelled by R K + K - and Rhad- For hadrons, real 
emission is restricted to pions only. 





KLOE 


BES 


BaBar 


&S+V, Rn+Ti- 


-1.36 


-0.818 


-0.0533 


&S+V, Rhad 


-1.06 


-1.81 


-0.1888 


OS+V+H, Rn+Ti- 


-0.186 


-0.0447 


-0.00229 


&S+V+H, Rhad 


0.47 


-0.15 


-0.0088 



IL(s) 
0.05 

0.04 

0.03 

0.02 

0.01 



-0.01 
-0.02 
-0.03 
-0.04 
-0.05 




o 2 

q 2 [GeV 2 } 



Figure 2.8: Hadronic vs. pion pair contributions to the real part of the vacuum polarisation function. 



the virtual and real contributions, so the careful evaluation of the integrals is needed. 
2.7 Exact NNLO numerical results 

In this Section we collect exact NNLO results which can serve as a benchmark for further investiga- 
tions. The event selections used here are very close to event selections used at meson factories to 
measure their luminosity and are described in Appendix A. We give separately the contributions 
from electron, muon, tau and pion pair production as well as the complete hadronic contribution. 
For the last one, as mentioned already in the previous Section, there exist no generator to give the 
contributions from the real hadron emission beyond the pion pair production. An educated guess 
of the size of the missing contributions, based on the fact that the pions are the lightest hadrons 
produced and that the highest energy of the meson factories is about 10 GeV, is that they should not 
be much higher than the contributions from pion pairs. Thus based on information from Table 2.7 



we can conclude that they should be completely negligible for the event selections used at meson 
factories for the luminosity measurements. 

In Tables 2.4-2.8 the meaning of the different entries is the following: a« is the Born cross 
section, a v+s the cross section with NNLO virtual plus soft photon corrections (Section 2.1), C/, 
the NNLO cross section with a self energy insertion corrected by the emission of one hard photon 
(Section 2.2), G v + s +h the sum of a v+s and a), and a pa i rs the leading-order cross section with emission 
of real pairs (Sections 2.3-2.5). The total NNLO massive correction can be obtained by summing 
Ov+s+h with Opairs ■ 

Table 2.4: Results for electron pair corrections at different energies, in GeV, for the reference 
event selection defined in Appendix A. The Ob is the Born cross section within the acceptance cuts. 
The separation of soft and hard photons is fixed at ffl/£beam = 10~ 4 , where ft) = £™ n (soft photon 
cut-off). All cross sections are given in nb. 



e + e 




Ob 


a h 


Gv+s 


&v+s+h 


Opairs 


KLOE 


1.020 


542.663(6) 


9.5021(2) 


-11.5666 


-2.0645(2) 


0.2712(15) 


BES 


3.097 


173.98(2) 


4.16202(13) 


-4.71708 


-0.55506(13) 


0.19977(116) 


BES 


3.650 


129.0958(4) 


3.19544(9) 


-3.55544 


-0.36000(9) 


0.188856(997) 


BES 


3.686 


123.32(1) 


3.14439(9) 


-3.49579 


-0.35140(9) 


0.18740(99) 


BaBar 


10.56 


5.481(1) 


0.202439(7) 


-0.223667 


-0.021228(7) 


0.01355(8) 


Belle 


10.58 


6.73555(4) 


0.21572(7) 


-0.25596 


-0.04024(7) 


0.0130999(469) 



Table 2.5: Results for muon pair corrections at different energies, in GeV, for the reference event 
selection defined in Appendix A. The Ob is the Born cross section within the acceptance cuts. The 
separation of soft and hard photons is fixed at ft)/.Ebeam = 10~ 4 , where ft) = E™ m (soft photon 
cut-off). All cross sections are given in nb. 







Ob 


O/, 


Oy +S 


Gv+s+h 


Opairs 


KLOE 


1.020 


542.663(6) 


1.49406(3) 


-1.7356(2) 


-0.2415(2) 


0.246(7)- 10" 7 


BES 


3.097 


173.98(2) 


1.01652(3) 


-1.09665(1) 


-0.08013(3) 


0.001337(5) 


BES 


3.650 


129.0958(4) 


0.83245(2) 


-0.88149(1) 


-0.04904(2) 


0.002003(6) 


BES 


3.686 


123.32(1) 


0.82215(2) 


-0.86988(1) 


-0.04773(2) 


0.002035(6) 


BaBar 


10.56 


5.481(1) 


0.075789(2) 


-0.079231(2) 


-0.003442(3) 


0.000451(2) 


Belle 


10.58 


6.73555(4) 


0.080377(8) 


-0.09009(1) 


-0.00971(1) 


0.0007587(14) 



As one can see from Tables 2.4-2.8 the discussed corrections cannot be ignored for high 
precision luminosity measurements. Actually the relative size of electron pair corrections amount 
to about 0.3% at KLOE, 0.1 - 0.2% at BES and BaBar, and below the 10 3 level at Belle only. 
The contribution of muon pair and hadronic corrections is slightly smaller, in the 0.05-0.1% range 
at all meson factories, while the contribution of tau pair corrections is, not surprisingly, generally 
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Table 2.6: Results for tau pair corrections at different energies, in GeV, for the reference event 
selection defined in Appendix A. The Ob is the Born cross section within the acceptance cuts. The 
separation of soft and hard photons is fixed at (o/E^ e3m = 10~ 4 , where (0 = E™ m (soft photon 
cut-off). All cross sections are given in nb. 







o B 


o h 


O v +s 


Ov+s+h 


Cpairs 


KLOE 


1.020 


542.663(6) 


0.0201637(4) 


-0.023412(2) 


-0.003248(2) 





BES 


3.097 


173.98(2) 


0.049672(2) 


-0.0540(1) 


-0.0044(1) 





BES 


3.650 


129.0958(4) 


0.058674(2) 


-0.0633(1) 


-0.0046(1) 





BES 


3.686 


123.32(1) 


0.057923(2) 


-0.0622(1) 


-0.0043(1) 





BaBar 


10.56 


5.481(1) 


0.0138398(4) 


-0.0144654(2) 


-0.0006257(5) 


0.120(3) -10- 8 


Belle 


10.58 


6.73555(4) 


0.014428(4) 


-0.01602(1) 


-0.00159(1) 


0.0000321(1) 



Table 2.7: Results for pion pair corrections at different energies, in GeV, for the reference event 
selection defined in Appendix A. The Ob is the Born cross section within the acceptance cuts. The 
separation of soft and hard photons is fixed at (o/E^ e3m = 10~ 4 , where ft) = E™ m (soft photon 
cut-off). All cross sections are given in nb. 



n + n 




o B 


Oh 


O v +s 


O v +s+h 


Op airs 


KLOE 


1.020 


542.663(6) 


1.17402(8) 


-1.35988(2) 


-0.18586(8) 





BES 


3.097 


173.98(2) 


0.95919(3) 


-1.03394(3) 


-0.07475 


0.000153(2) 


BES 


3.650 


129.0958(4) 


0.77337(2) 


-0.81806(3) 


-0.04469 


0.000539(7) 


BES 


3.686 


123.32(1) 


0.76286(2) 


-0.80626(3) 


-0.04340 


0.000564(8) 


BaBar 


10.56 


5.481(1) 


0.051037(2) 


-0.053328(3) 


-0.002291(4) 


0.000029(3) 


Belle 


10.58 


6.73555(4) 


0.054457(6) 


-0.0612(1) 


-0.0067(1) 


0.00015(1) 



negligible.Therefore, the pair corrections have to be implemented into a MC event generator at 
least in an approximated form (see the next Section for their implementation in the generator 
BabaYaga@NLO). In particular, for the real emission one can conclude that only the reaction 
e + e~ — > e + e~e + e~ gives significant contributions to the cross section used in the luminosity 
measurements. When the accuracy of the experiment reaches the level 10~ 3 this process has to be 
considered and its contributions added to the theoretical cross section or alternatively subtracted as a 
background from the experimental cross section. 

3. The NNLO massive corrections in BabaYaga@NLO 
3.1 The program 

BabaYaga is an event generator for precise simulations of the processes e + e~ — > e + e~ , , yy in 
QED. It was developed for precision measurements with per mille accuracy of the luminosity of 
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Table 2.8: Results for hadronic corrections at different energies, in GeV, for the reference event 
selection defined in Appendix A. The Ob is the Born cross section within the acceptance cuts. The 
separation of soft and hard photons is fixed at (o/E^ cam = 10~ 4 , where (0 = E™ m (soft photon 
cut-off). All cross sections are given in nb. 



hadrons 




o B 


Oh 


Ov+s 


Oy+s+h 


Cpion pair only 


KLOE 


1.020 


542.663(6) 


1.5248(6) 


-1.062(8) 


0.463(8) 





BES 


3.650 


129.0958(4) 


1.66065(8) 


-1.81(1) 


-0.15(1) 


0.000539(7) 


BaBar 


10.56 


5.481(1) 


0.17995(2) 


-0.1888(4) 


-0.0088(4) 


0.000029(3) 


Belle 


10.58 


6.73555(4) 


0.18969(1) 


-0.2124(5) 


-0.0227(5) 


0.00015(1) 



GeV scale e + e colliders. It has been adopted and is still presently used for this purpose by KLOE, 
BES/BES-III, CLEO, BaBar and Belle experiments. 

BabaYaga is released in two versions, the more accurate BabaYaga@NLO and BabaYaga 
3.5, http://www2.pv.infn.it/~hepcomplex/babayaga.html. In BabaYaga 3.5 the most 
relevant QED radiative corrections are included in a pure Parton Shower (PS) approach, even though 
improved to include radiation interference effects for a more accurate description of radiative events. 
BabaYaga @NLO also includes non-logarithmically enhanced O(a) corrections, which were the 
main source of theoretical error of BabaYaga 3.5. As summarised in Section 3.2, the necessary 
NLO ingredients are matched in BabaYaga @NLO with the PS approach, in order to preserve 
exponentiation of large contributions and ensure normalisation at NLO accuracy. 

Concerning the NNLO corrections that are the concern of the present study, the contribution of 
NNLO massive corrections included in the code is approximate and comes, as detailed in Section 
3.3, from 

a) insertion of self-energy corrections in NLO virtual + soft photon correction; 

b) insertion of self-energy corrections in NLO hard photon correction. 

From a theoretical point of view, the part b) coincides with the contribution to the NNLO 
corrections described in Section 2.2, although the e + e~ — > e + e~y matrix element, phase space and 
related MC integration are completely independent of the calculation previously discussed. From 
a numerical point of view, one should expect a priori that the two calculations provide results in 
agreement within the respective statistical uncertainties, as it will be studied in the following. On 
the other hand, the contribution a) is just a subset of the full NNLO correction described in Section 
2.1, as further detailed in Section 3.3. 

Furthermore, the contributions due to real pair emission of the type 2 — > 2 + (2), which are 
included in the complete NNLO benchmark calculation (Sections 2.3-2.5), are presently neglected 
in the program. 

3.2 General formulation of BabaYaga @ NLO 

As far as photon corrections are concerned, the two basic ingredients combined in BabaYaga® NLO 
to guarantee the target theoretical accuracy are 
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1. exact NLO QED corrections (soft+ virtual and hard contributions); 



2. leading QED logarithms due to multiple collinear and soft radiation beyond 0(a). 
The matching is performed according to the following formula [31] 

°° 1 ( " \ 

^matched = F SV n(Q 2 ,S) £ - JJFhj \^ LL \ 2 d<$> n (3.1) 



n=0 n - \i=0 



where 



1. n(<2 2 ,£) is the Sudakov form factor. It describes universal (process independent) virtual + 
soft radiation up to the energy fraction e (soft-hard separator) 

U(Q 2 ,e) = exp ( -— /+ If) , JJ = log^r, /+ = / dzP(z) 
P(z) 



2tc 1 / ' m? 

1+z 2 
l-z ' 







In the program, the scale Q 2 is chosen such that L' = log % = log — 1 = L — 1 where s, t 
and u are the Mandelstam variables of the process and m e is the electron mass. This choice 
is dictated by the perturbative NLO calculation of the Bhabha process and ensures that, in 
addition to initial and final state leading contributions, also initial-final state interference 
leading effects are resummed to all orders. 

2- \^n,u\ • «-photon emission squared amplitude in collinear approximation, with phase space 
factor <i<J>„ 

3. Fsv and Fu'- residuals of the exact NLO calculation w.r.t. the leading log approximation 
ensured by the ingredients 1. and 2. above, i.e. 

l| 2 -|^i,ii| 2 



Fsv — 1 + (C a — C a LL) , Fh — H j— -7. r~ 

C a : exactsoftplus virtualNLO^factor, 

C a ,LL '■ 0(a) expansion of the Sudakov form factor 

i| 2 : exact NLO hard bremsstrahlung squared matrix element 



The explicit expressions of all the above contributions can be found in [31]. 

A further necessary ingredient is the correction due to vacuum polarisation. It is included in the 
Bhabha Born matrix element setting r s = a(s)/a and r t = a(t) / a and rescaling the s and t channel 
amplitudes as follows 

|^b| 2 = \-^Oj + -^oA 2 -> l-Av/>| 2 = \^o, s r s + ^o, t r t \ 2 ■ ( 3 - 2 ) 

In the code, we use the resummed expression a(q 2 ) = a/(l —Aa(q 2 )), where Aa(q 2 ) is the 
fermionic contribution to the photon self-energy. It is treated analytically for the leptonic and 
top-quark one-loop contributions, while the non-perturbative five quark (hadronic) contribution, 
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AoCi wd , is included according to the latest Jegerlehner [50] and Teubner et al. [47] parameterisations. 
In order to include an important class of 0(a 2 ) factorizable corrections, we insert the vacuum 
polarisation correction in the NLO cross section too, both in the soft plus virtual contribution and the 
hard photon matrix element. The inclusion of the vacuum polarisation both in the soft plus virtual 
and hard photon part guarantees that their sum is independent of the soft-hard separator e, as we 
explicitly checked. 



3.3 NNLO massive corrections in BabaYaga@NLO 

As detailed in [31], it is possible to extract from the matched formula given in Eq. (3.1) the different 
pieces contributing to the cross section at NNLO. In particular, to explain how NNLO massive 
corrections are (approximately) taken into account in BabaYaga@NLO, let us write the expansion 
up to 0(a 2 ) of the cross section with soft plus virtual corrections. It can be derived from the 
first (n = 0) term of the infinite sum in Eq. (3.1). In order to highlight the s, t and interference 
contributions, first we define 



do = do S fi + do,fl + do st fi 
= (B s + B t + B st )do , 
dog v = do? sv + do? sv + doZ sv 
= (E s + E t + E st )da Q , 



(3.3) 



where B s , B t and B st are the percentage s, t and st lowest-order contributions to the complete 
Bhabha differential cross section <f Go, and E s , E t and E st are the NLO SV correction factors for 
each contribution, in units of dOQ. Truncating every factor in Eq. (3.3), improved with vacuum 
polarisation effects as described above, we get from Eq. (3.1) at 0(a 2 ) 



dosv 
d<Jo 



1+V + 



[1 + (E s - VB s )r 2 + (E t - VB t )rf + (E st - VB st )r s r t ] 
{B s r 2 s +B t r 2 +B st r s r t ), 



(3.4) 



where V = —(2a/n)I + L' is the 0(a) truncation of the Sudakov form factor, Ej and B,- have been 
defined above and r Sjt are the vacuum polarisation corrections for the s and t channels, respectively. 
If we define 1/(1 - /S.a(q 2 )) = 1/(1 — S q i), the r§, rf and r s r t read 

r 2 = 1 +25, + 35/, 
r 2 = l+25 f + 35, 2 , 
r s r t = \ + 8 s + 8 t + S 2 + 8 2 + S s S t . (3.5) 
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Retaining only terms up to 0(a 2 ), Eq. (3.4) reads 

dosv _ l 
do 



+ V + (E s - VB S ) + (E t - VB t ) + (Est ~ VB st ) 
+ 2(B.,8 s +B t 8 t )+B st (8 s + 8 t ) 
+ 1/2V 2 

+ (E s - VB S )8 S + (E t - VB t )8 t + (E st - VB S ,)(8 S + 8 t ) 

+ 3(B S 8 S 2 +B t 8 2 ) +B a (Sf + 8? + 8 s 8 t ) 

+ V [(E s - VB S ) + (E t -VB t ) + (E st - VB st ) 

+ V[2(B s 8 s +B t 8 t )+B st (8 s + 8,)} 

+ [(E s - VB S ) + (E t - VB t ) + (E st - VB^)} 

x [2(B s 8s+BtSt)+Bs t (S s + S t )]. 



(3.6) 



The first line of Eq. (3.6) is the Born contribution, the second line is the one soft photon plus 
one loop virtual correction (notice that it is equal to E s + E t +E st because B s + B t +B st = 1), the 
third line is the vacuum polarisation correction at 0(a) and the remaining lines represent the cross 
section with 0(a 2 ) soft plus virtual corrections. From the latter it is simple to disentangle the 
NNLO contribution due to the insertion of the vacuum polarisation correction in the 0(a) soft plus 
virtual coefficients. The relevant terms are those containing a 5; factor. Among them, there is a 
pure two-loop self-energy correction (sixth line in Eq. (3.6)), which we discard for the comparison 
with the exact NNLO calculation, in accordance with the discussion of Section 2. 1 . Therefore the 
formula of interest reduces to 



Equation (3.7) is used in the present study to validate the approximate treatment of NNLO massive 
corrections as in BabaYaga@NLO in the soft plus virtual regime. 

Equation (3.7) must be added to the contribution obtained by dressing the hard bremsstrahlung 
cross section with self-energy corrections. The hard photon matrix element is the sum of eight 
amplitudes where the real photon is attached to a s or t channel-like diagram. As in Eq. (3.2), those 
amplitudes are rescaled by r s and r t , respectively, to account for the effect of vacuum polarisation. 
In BabaYaga@NLO, the squared amplitude for the emission of a real photon \^\\ is exact as 
calculated through FORM [51] and cross checked with the output of the ALPHA algorithm [52]. 

To summarize, the best knowledge of running a QED and all the virtual factorizable NNLO 
vacuum polarization corrections are included in BabaYaga@NLO. In particular, within the full 
set of (reducible and irreducible) NNLO virtual massive corrections present in the exact calculation 
described in Section 2.1, only the subset of loop-by-loop corrections (see Figure 2.2) is taken into 
account in the code. In other words, the contributions of purely irreducible NNLO vertex and box 



dosv 
d(jQ 



(E s - VB s )Ss + (E t - VB t )8 t + (E^ - VB st )(8 s + 8 t ) 

+ V[2(B s 8 s +B t 8 t )+B st (8 s + 8,)] 

+ [(E s - VB S ) + (E t - VB t ) + (E^ - VB^] 

x [2(BsSs+BtSt)+Bst(Ss + 8t)]. 



(3.7) 



corrections is not included. However, since the real pair emission corrections are neglected as well, 
this means that there is no imbalance of ln 3 (s/m e ) terms in BabaYaga@NLO, thus respecting the 
compensation mechanism of leading mass singularities discussed in Section 2.1 and Section 2.3 
(see Eq. (2.27) and Eq. (2.37)). 

3.4 Numerical results of BabaYaga@NLO 

Below we give the benchmark results from the BabaYaga@NLO MC event generator. Their 
detailed comparison to the exact results will be given in the next Section. 

Table 3.1: B ABAYAGA @ NLO results for electron pair corrections at different energies, in GeV, 
for the reference event selection defined in Appendix A. The Ob is the Born cross section within 
the acceptance cuts. The separation of soft and hard photons is fixed at (o/E\, eam = 10~ 4 , where 
CO = E™ m ( soft photon cut-off). All cross sections are given in nb. 



e^e 




o B 


Oh 


O v+s 


Ov+s+h fpairs 


KLOE 


1.020 


542.663(6) 


9.5022(8) 


-11.0721(4) 


-1.5699(9) 


BES 


3.097 


173.98(2) 


4.1624(4) 


-4.4818(2) 


-0.3194(5) 


BES 


3.650 


129.0958(4) 


3.1960(3) 


-3.3730(2) 


-0.1770(4) 


BES 


3.686 


123.32(1) 


3.1447(3) 


-3.3163(2) 


-0.1716(4) 


BaBar 


10.56 


5.481(1) 


0.20244(2) 


-0.20971(5) 


-0.00727(5) - 


Belle 


10.58 


6.73555(4) 


0.21563(2) 


-0.23994(2) 


-0.02431(3) - 



Table 3.2: BABAYAGA@NLO results for muon pair corrections at different energies, in GeV, 
for the reference event selection defined in Appendix A. The Ob is the Born cross section within 
the acceptance cuts. The separation of soft and hard photons is fixed at w/isbeam = 10~ 4 , where 
CO = E^ m ( soft photon cut-off). All cross sections are given in nb. 



MM 




o B 


Oh 


o v+s 


Ov+s+h ^pairs 


KLOE 


1.020 


542.663(6) 


1.4942(2) 


-1.7441(2) 


-0.2499(3) 


BES 


3.097 


173.98(2) 


1.01672(9) 


-1.0960(2) 


-0.0793(2) 


BES 


3.650 


129.0958(4) 


0.83252(7) 


-0.88041(9) 


-0.0479(1) 


BES 


3.686 


123.32(1) 


0.82221(7) 


-0.8688(1) 


-0.0466(1) 


BaBar 


10.56 


5.481(1) 


0.07580(1) 


-0.07872(2) 


-0.00292(2) - 


Belle 


10.58 


6.73555(4) 


0.080376(6) 


-0.08948(2) 


-0.00910(2) - 



We repeat here that BabaYaga@NLO does not generate the real pair emission contribution. 
This explains why the corresponding predictions for the cross section a pa i rs and Ohadrons are not given 
in Tables 3.1-3.4. However, from Section 2.7 (see the next Section for a more detailed discussion) it 
is clear that only the reaction e + e~ — > e + e~e + e~ gives a contribution which is relevant and needs a 



-24- 



Table 3.3: BABAYAGA@NLO results for tau pair corrections at different energies, in GeV, for 
the reference event selection defined in Appendix A. The Ob is the Born cross section within the 
acceptance cuts. The separation of soft and hard photons is fixed at CO/E\, eam = 10~ 4 , where 
CO = E™ m ( soft photon cut-off). All cross sections are given in nb. 



T T y/s Ob Oj t O v+s Ov+s+h Cpairs 

KLOE 1.020 542.663(6) 0.020166(3) -0.023704(2) -0.003538(4) 

BES 3.097 173.98(2) 0.049683(5) -0.05421(1) -0.00453(1) 

BES 3.650 129.0958(4) 0.058679(7) -0.06323(2) -0.00455(2) 

BES 3.686 123.32(1) 0.057928(7) -0.06219(2) -0.00426(2) 

BaBar 10.56 5.481(1) 0.013847(4) -0.014541(4) -0.000694(6) - 

Belle 10.58 6.73555(4) 0.014423(1) -0.016091(7) -0.001668(7) - 

Table 3.4: BABAYAGA@NLO results for hadronic corrections at different energies, in GeV, for 
the reference event selection defined in Appendix A. The Ob is the Born cross section within the 
acceptance cuts. The separation of soft and hard photons is fixed at CO/E\, eam = 10~ 4 , where 
CO = E^ m ( soft photon cut-off). All cross sections are given in nb. 



hadrons y/s Og Oh OV-M Gy+s+h Ohadrons 

KLOE 1.020 542.663(6) 1.5247(5) -1.126(2) 0.399(2) 

BES 3.650 129.0958(4) 1.6613(3) -1.7860(2) -0.1247(4) 

BaBar 10.56 5.481(1) 0.17984(2) -0.18760(4) -0.00776(5) 

Belle 10.58 6.73555(4) 0.18964(3) -0.21089(5) -0.02125(6) 



particularly careful treatment. To account for this reaction a dedicated generator should be used (for 
example Helac-Phegas). 

Comparing the exact NNLO results of Tables 2.4-2.8 with those of BABAYAGA@NLO one 
can note that, while differences are present as expected, in the soft + virtual cross section, the 
independent predictions of the two calculations for the hard photonic correction a/, are in agreement 
within the MC statistical uncertainty. 

4. BabaYaga@NLO versus the exact NNLO massive corrections 

In this Section we would like to answer how well the NNLO massive corrections are accounted for in 
the BabaYaga@NLO event generator. The results are summarised in Table 4.1 and Figs. 4.1-4.4. 

In Table 4.1 we show the impact of the exact corrections and the corrections given by 
BabaYaga@NLO for the reference event selections at meson factories. One can observe that the 
leptonic corrections are dominated by the electron corrections. Moreover they are well accounted 
for in the BabaYaga@NLO code with the biggest unaccounted correction of the order of 0.5 %c. 
In Figs. 4.1-4.4 we have studied the stability of the obtained results against changes of the event 
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Table 4.1: Comparison of the exact massive NNLO with BABAYAGA@NLO results given for 
the reference event selections at different meson factories. The event selections are defined in the 
Appendix A. The a BY is the cross section in nb from BabaYaga@NLO, and S x = -5— with 
x = e + e~ ,lep,tot, where tot stands for leptonic (lep) + hadronic corrections. 









Oby 


Se+e- 


Slep [%o] 


Shad [%°] 


Stot ]%o\ 


KLOE 


1.020 


NNLO 




-3.935(4) 


-4.472(4) 


1.02(2) 


-3.45(2) 






BabaYaga@NLO 


455.71 


-3.445(2) 


-4.001(2) 


0.876(5) 


-3.126(5) 


BES 


3.097 


NNLO 




-2.246(8) 


-2.771(8) 










BabaYaga@NLO 


158.23 


-2.019(3) 


-2.548(3) 






BES 


3.650 


NNLO 




-1.469(9) 


-1.913(9) 


-1.3(1) 


-3.2(1) 






BabaYaga@NLO 


116.41 


-1.521(4) 


-1.971(4) 


-1.071(4) 


-3.042(5) 


BES 


3.686 


NNLO 




-1.435(8) 


-1.873(8) 










BabaYaga@NLO 


114.27 


-1.502(4) 


-1.947(4) 






BaBar 


10.56 


NNLO 




-1.48(2) 


-2.17(2) 


-1.69(8) 


-3.86(8) 






BabaYaga@NLO 


5.195 


-1.40(1) 


-2.09(1) 


-1.49(1) 


-3.58(2) 


Belle 


10.58 


NNLO 




-4.93(2) 


-6.84(2) 


-4.1(1) 


-10.9(1) 






BABAYAGA@NLO 


5.501 


-4.42(1) 


-6.38(1) 


-3.86(1) 


-10.24(2) 



selection. We plot there the relative difference between the exact massive NNLO corrections and the 

jj.NNLO_g.NNLO 

BabaYaga@NLO results, i.e. '"" gBy BY — , where a BY is the full (with all radiative corrections 
included) prediction of BabaYaga@NLO according to Eq. (3.1). One can observe that even if 
for very stringent cuts one can have the differences up to 0.9 %c (Belle, Fig.4.3), generally the 
results are not varying rapidly. Comparisons between results for Belle (Fig.4.3) and BaBar event 
selections (Fig.4.4) show also that the actual difference is sensitive to the event selection used and it 
is recommended to study the effect on the NNLO massive corrections, whenever the event selection 
is changed by a given experiment. Similar effects are observed for the hadronic corrections (with 
the exception of J/y and y(2S) energies, which will be examined in a separate paper). The biggest 
hadronic correction missing in BabaYaga@NLO is about 0.4 per mille and in addition for KLOE 
and BES energies the missing hadronic contribution is of the opposite sign of the missing leptonic 
contribution, thus partly cancelling each other. For B-factories the missing leptonic and hadronic 
corrections are of the same sign, but even there the sum of the missing parts does not exceed one per 
mille (0.7 per mille for the reference event selections). 

5. Conclusions 

We presented an exact calculation of NNLO massive corrections to Bhabha scattering, beyond ap- 
proximations previously addressed in the literature, and performed detailed studies of the missing part 
of these corrections in the BabaYaga@NLO MC program. The approximate BabaYaga@NLO 
formulae were confronted with the exact NNLO results for event selections close to the experimental 
ones. The stability of the results with changing of the event selections was also examined. 
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Figure 4.1: T/ie relative difference gBY BY ofNNLO massive leptonic and hadronic corrections 
between exact and BABAYAGA@NLO, as a function of acollinearity cut for two different angular 
acceptance regions for KLOE-like event selections (see Appendix A). 



We found that NNLO massive conections are relevant for precision luminosity measurements 
with 10~ 3 accuracy, their relative total contribution to the Bhabha scattering cross section being 
of a few per mille. Hierarchically, the electron pair contribution is the largely dominant part of 
the correction, the muon pair and hadronic conection, which are the next-to-important effect, are 
quantitatively on the same grounds, while the tau pair contribution is negligible for the energies of 
meson factories. 

Thanks to the exact NNLO calculation we tested the theoretical accuracy of the generator 
BabaYaga@NLO which includes such conections according to an approximate treatment through 
insertion of self-energy contributions into NLO conection factors. On the grounds of several 
numerical results, we concluded that the very bulk of the conection is taken into account in the 
program. For reference realistic event selections the maximum observed difference is at the level 
of 0.07%. When cuts are varied the sum of the missing pieces can reach 0.1%, but for very tight 
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Figure 4.2: 77ie relative difference ^ gBY BY ofNNLO massive leptonic and hadronic corrections 
between exact and BABAYAGA@NLO, as a function of an angular cut for BES-like event selections 
( see Appendix A ). 



acollinearity cuts only. 

As a possible perspective, it is worth mentioning that the leading logarithmic part of the missing 
pieces in BabaYaga@NLO, coming from the interplay between NNLO virtual and real pair 
corrections, could be accounted for by means of appropriate singlet/non-singlet QED structure 
functions, as discussed e.g. in [53, 54, 55] and done in the past for small-angle Bhabha scattering 
at LEP [56]. However, this improvement of the theoretical formulation of BabaYaga@NLO 
would be necessary only whenever the experimental precision of the luminosity measurements 
should require it and should be accompanied by the inclusion of other sources of NNLO ingredients 
presently neglected in the code and contributing at an accuracy below the 10~ 3 level [1, 57]. As 
already discussed in the paper, open issues of the present work are a more detailed study of hadronic 
NNLO corrections and of the uncertainty induced by hadronic vacuum polarisation insertions to 
Bhabha scattering in a close vicinity of nanow resonances. We plan to address the above perspectives 
in future works. 
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A. The experimental cuts 

A.l <I> factory KLOE/DA<I>NE (Frascati) - the reference event selection 

• y/s= 1.02 GeV 

• £ min = 0.4 GeV 

• 55° < 0± < 125° 

• the maximal allowed 2D acollinearity of two charged tracks <^ max = 9° 

To see how the NNLO corrections depend on the event selection we obtained results also for 

i. wider selection 20° < 0± < 160° 

ii. For each 0± selection the following £ max was used (£ max = 4° , 5° , 6° , 7° , 8° , 9° , 10° , 1 1° , 12° , 13° , 14°) 

A.2 charm/r factory BES-III (Beijing) - the reference event selection 

• y/s = 3.686 GeV, 3.65 GeV and 3.097 GeV 

• | cos 0| < 0.8, where is the polar angle of the electron or positron in the lab system, this 
corresponds to the barrel region of BES-III detector. Since in BEPC, e + and e~ beams are 
colliding with equal energy but at a 22mrad crossing angle, the lab system is slightly different 
from the CoM system. 

• E e + > 1.0 GeV and E e - > 1.0 GeV, where E is the energy deposited in the electromagnetic 
calorimeter (EMC). 1 

To see how the NNLO corrections depend on the event selection we obtained results also for 
| cos 1 <0.7, 0.75, 0.85 and 0.9. 

A.3 B factory BaBar/PEP-II (SLAC) - the reference event selection 

• y/s = 10.56 GeV 

• \p + \/E beam > 0.75 and \p-\/E bsam > 0.50 
or \p-\/E beam > 0.75 and \p+\/E heam > 0.50 

• |cos(0±)| <0.65 and |cos(0+)| < 0.60 or |cos(0-)| < 0.60 

• the maximal allowed 3D acollinearity of two charged tracks <^3dmax = 30° 

To see how the NNLO corrections depend on the event selection we obtained results also for 

i. |cos(0±)| < 0.70 and |cos(0+)| < 0.65 or |cos(0-)| < 0.65 

ii. |cos(0±)| < 0.60 and |cos(0+)| < 0.55 or|cos(0-)| < 0.55 
For each | cos(0±)| selection we have used the following <^3dmax 
(fedmax = 20° , 22° , 24° , 26° , 28° , 30° , 32° , 34° , 36° , 38° , 40° ) 



e + and e deposit virtually all their energy in the EMC, so this is the energy they carried. 
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A.4 B factory Belle (KEK) - the reference event selection 

Belle runs at an asymmetric e + e~ collider, but all criteria are expressed in the CoM frame. The 
selection of Bhabha events by Belle are: 

• ^=!0-58GeV 

• 50.5° < 6± < (180-50.5)° 

• Two charged tracks have momentum > 2.645 GeV 

• The track with maximum deposited energy in EMC greater than 2 GeV, 

• The sum of the deposited energies of all tracks in EMC is greater than 4 GeV (both charged 
and neutral particle can deposit energy in EMC and it is not checked if the particle is charged 
or neutral) 

• Acollinearity angle (2D) <^2dmax < 10° 

• Transverse momentum of any observed charged particle greater than 0. 1 GeV 

To see how the NNLO corrections depend on the event selection we obtained results also for 

i. 45.5° < 6± < (180-45.5)° 

ii. 55.5° < 0± < (180-55.5)° 

For each G± selection we have used the following § 2 dmax (fedmax = 5°, 6°, 7°, 8°, 9°, 10°, 11°, 12°, 13°, 14°, 15°) 
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